Dynamical Mean Field Theory of Quantum Stripe Glasses 
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We present a many body approach for non-equilibrium behavior and self-generated glassiness in 
strongly correlated quantum systems. It combines the dynamical mean field theory of equilibrium 
systems with the replica theory for classical glasses without quenched disorder. We apply this 
approach to study a quantized version of the Brazovskii model and find a self-generated quantum 
glass that remains in a quantum mechanically mixed state as T — > 0. This quantum glass is formed 
by a large number of competing states spread over an energy region which is determined within our 
theory. 

PACS numbers: 61.43. Gt, 75.10.Nr, 74.25.2q 



S3 
S3 



I 

S3 
O 

o 



> ■ 
a\ ■ 
o\ : 
o . 

(N ■ 

o ■ 
<n : 
o . 

■ 



I 

S3 
O 

o 



X 



I. INTRODUCTION 

The formation of glasses upon cooling is a well known 
phenomenon for classical liquids. Even without quenched 
disorder, the relaxation times become so large that a 
frozen non-ergodic state is reached before the nucleation 
into the crystalline solid sets in. The nucleation is es- 
pecially easy to avoid if the material has many poly- 
morphisms as does for example SiCb- The system be- 
comes unable to reach an equilibrium configuration on 
the laboratory time scale and exhibits aging and mem- 
ory effects. While extrinsic disorder is widely acknowl- 
edged to lead to glassy phenomena in quantum systems, 
can a similar self generated glassiness occur for quan- 
tum liquids? Candidate materials for this behavior are 
strongly correlated electron systems which often exhibit 
a competition between numerous locally ordered states 
with comparable energies. Examples for such behav- 
ior are colossal magneto- resistance materials ; 1 ^ 3 cuprate 
superconductors^iSiLSiSiiSiil and likely low density elec- 
tron systems^ possibly close to Wigner crystallization. 

The possibility of self generated quantum glassiness 
in such systems and the nature of the slow quantum 
dynamics has been little explored. Within a purely 
classical theory a "stripe glass" state was recently pro- 
posed in Refs i 13 i 14 . This proposal was stimulated by 
nuclear magnetic resonance^S*^ and /it-spin relaxation 
experiment oSi 10 ! 11 who found static or quasi-static charge 
and spin configurations similar to glassy or disordered 
systems. A prominent effect which reflects this slowing 
down was the "wipe-out" effect where below a certain 
temperature the typical time scales become so long that 
the resulting rapid spin lattice relaxation cannot be ob- 
served anymor o 5 i 6 i 8 . More recently, the dynamical be- 
havior of the stripe glass of RefsiA2i±i, as determined by 
Grousson et alw^ was found to be in quantitative agree- 
ment with NMR-experimentsi^. Finally, recent /z-spin 



experiment: 
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analyzed the freezing temperatures as 



function of charge carrier concentration and disorder con- 
centration and found a quantum glass transition which 
was insensitive to the amount of disorder added to the 



system. This latter experiment, which seems to support 
a generic explanation for glassiness, not caused by impu- 
rities, also demonstrated the need for a more detailed 
investigation of the quantum regime of glassy systems. 
It is then important to develop appropriate theoretical 
tools to predict whether a given theoretical model for a 
quantum many body system will exhibit self-generated 
glassiness where the system forms a glass for arbitrary 
weak disorder. 

In this paper we develop a general approach to self 
generated quantum glasses that combines the dynami- 
cal mean field theory (DMFT) of quantum many body 
systemsiiiSiifi with the replica technique of Refs l2fl2ll 
which describes classical glasses without quenched dis- 
order. In addition, we also go beyond strict mean field 
theory and estimate the spectrum of competing quantum 
states after very long times. This approach is applied 
to study a model with competing interactions which has 
a fluctuation induced first-order transition to a striped 
(2D) / lamellar (3D) phase. A quantum glass is shown to 
be in a quantum mechanical mixed state even at T = 0, 
formed by a large number of states which can be consid- 
erably above the true ground state. Using the concept 
of an effective temperatur e 22 ^^ 2i for the distribution of 
competing ground states, our theory allows the investi- 
gation of glassy non-equilibrium dynamics in quantum 
systems using standard techniques of equilibrium quan- 
tum many body theory. On the mean field level, our 
theory is in complete agreement with the explicit dynam- 
ical non-equilibrium approach developed by Cugliandolo 
and Lozano 2 ^ for quenched disordered spin glasses with 
entropy crisis, which includes weak long-term memory 
effects and the subtle interplay of aging and stationary 
dynamics. 

We find that when quantum fluctuations are weak, i.e., 
for a small quantum parameter, c (for a definition see 
EqH below) , the glass transition resembles that of clas- 
sical models^ that exhibit a dynamical transition at a 
temperature Ta (where mean field theory becomes non- 
ergodic) and a Kauzmann entropy crisis at Tk < Ta if 
equilibrium were to be achieved. The actual laboratory 
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glass transition is located between these two tempera- 
tures and depends for example on the cooling rate of the 
system. Beyond a critical c, T K and Ta merge and the 
transitions change character. Within mean field theory, 
a discontinuous change of the relevant quantum mechan- 
ical states occurs at the glass transition. Even going be- 
yond mean field, the system remains in a mixed quantum 
state however with a non-extensive number of relevant 
quantum states. Still, at the quantum glass transition 
a discontinuous change of the density matrix to a usual 
quantum liquid occurs. Put another way, a glass in con- 
tact with a bath at T = is essentially a classical object, 
qualitatively distinct from a quantum fluid, enforcing the 
quantum glass transition to be discontinuous. Our results 
clearly support the later scenario as can be seen in Fig[21 
below. 

In classical glass forming liquids an excess or config- 
urational entropy with respect to the solid state due to 
an exponentially large number of metastable configura- 
tions emerges below Ta . Obviously, in the quantum limit 
a glass transition must be qualitatively different. Even 
a very large number of long lived excited states cannot 
compensate their vanishing Boltzmann weight at equi- 
librium for T = 0. An exponentially large ground state 
degeneracy on the other hand is typically lifted by hy- 
bridization, saving kinetic energy. In Ref|23 it was then 
argued that for a glassy quantum system the Edwards- 
Anderson order parameter vanishes continuously at the 
quantum glass transition (in contrast to the classical be- 
havior). On the other hand, in Refl2Sl it was concluded 
that in certain spin glasses with quenched disorder the 
glass transition at sufficiently small T is of first order. 

The outline of this paper is as follows. In section II we 
introduce a quantized version of the Brazovskii model of 
micro-phase separation. We then discuss the replica ap- 
proach used in this paper as well as the dynamical mean 
field theory employed for the solution of the problem in 
sections III and IV, respectively. Conclusions which are 
based on the approach developed here but which go be- 
yond the strict mean field limit are discussed in section 
V. Finally we present a summary of our results in the 
concluding section VI. The equivalence of the "cloned- 
liquid" replica approach and the Schwinger-Keldysh the- 
ory of non-equilibrium quantum systems is presented in 
the appendix. 



II. THE MODEL 

We consider a Bose system with field, <p x , governed 
by the action S [ip x ] with competing interactions which 
cause a glass transition in the classical limit. Specifically 
we consider a system with action 

+2, -»([V 2 + , 2 ]^) 2 ). 



Here, qo is a wave number which supports strong fluctu- 
ations for momenta with amplitude |q| = qo, i.e. mod- 
ulated field configurations. The classical version of the 
model, EqQ] was shown by Brazovskii-- to give rise to a 
fluctuation induced first order transition to a lamellar or 
smectic state. Within equilibrium statistical mechanics, 
this ordered state gives the lowest known free energy. In 
RefsOIl we demonstrated that, within non-equilibrium 
classical statistical mechanics, an alternative scenario is a 
self generated glass. Instead of the transition to a smec- 
tic state, metastable solutions built by a superposition of 
large amplitude waves of wavenumber go > but with ran- 
dom orientations and phases eme rge^. Those form the 
stripe glass state discussed in Ref. Il3ll4 It is unclear so 
far whether the glassy solution occurs only if the ordered 
phase can be avoided by super-cooling or whether there is 
a parameter regime where the glass is favored regardless 
of the cooling rate. 

In EqQ]we consider additional quantum fluctuations, 
characterized by the velocity, c. Clearly, quantum fluc- 
tuations reduce the tendency towards a fluctuation in- 
duced first order transition. This can be seen by eval- 
uating the (<y9 2 (a;,r)) within the spherical approxima- 
tion. In the classical limit (ip 2 (x, t)) ~ Tq^r^ 1 ^ 2 with 
renormalized mass r = r$ + u(ip 2 (x, r)) and the solu- 
tion r = clearly does not exist. The same fluctuations 
which suppress the occurrence of a second order transi- 
tion lead to a first order transition at the temperature 
where ro ~ uTq^r 1 / 2 ^£ In the quantum limit the be- 
havior is conceptually similar but fluctuations grow only 
logarithmically, (</? 2 (x,t)) ~ g 2 log(^). For exponen- 
tially large correlation length r -1 / 2 there should also be 
a fluctuation induced first order transition to a smectic, 
which might be related to the state proposed in Ref 1.3 ll in 
the context of strongly correlated quantum systems. An- 
other option however is the emergence of a stripe glass, 
even for large quantum fluctuations, which results in an 
amorphous modulated state instead. The investigation 
of this option will be the specific application of our the- 
ory. Before we go into specifics of the model, EqQJ we 
develop a general framework for the description of self 
generated quantum glasses. 



III. THE "CLONED LIQUID" - REPLICA 
APPROACH 

Competing interactions of a glassy system cause the 
ground state energy as well as the excitations to be very 
sensitive to small additional perturbations. In order to 
quantify this we introduce, following Refl2fl. a static "er- 
godicity breaking" field ip: 

9 f 3 2 

S^, [if] = S [tp] + - / drd x (V> x - </?x (t)) 

and take the limit 5^0 eventually. The coupling be- 
tween f and tp will bias the original energy landscape 
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in the "direction" of the configuration ip, enabling us 
to count distinct configurations. Adopting a mean-field 
strategy, we assume that even in the quantum limit ip 
should be chosen as static variable, probing only time 
averaged configurations. 

Introducing = — TlogZ^, with the biased partition 
function Z$ — J Dipe~ s ^ , f^{T — > 0) corresponds to 
the ground state energy for a given ip. If there are many 
competing ground states, it is natural to assume that 
determines the probability, p^, for a given configuration. 
If we identify the actual state of the system we gain the 
information S c = — lim s ^o f Dip p^ log p^. Maximizing 
this configurational entropy S c with respect to p^ yields 
the usual result 

oc exp (-/^/T e ff) , (2) 

where the effective temperature, T eff is the Lagrange mul- 
tiplier enforcing the constraint that the typical energy is 
F = lim s ^o / Dipp^fy. T e ff is a measure for the width of 
the energy region within which the relevant ground state 
energies can be found. If the system is glassy, T e g > T 
and a quantum mechanically mixed state results even as 
T — > 0. Thus, a glass will not be in a pure quantum me- 
chanical state (characterized by a single wave function) 
even at zero temperature. 

Introducing the ratio to = ^ it follows p^ oc Z™, 
leading to 

~ d (mF (to)) 
dm 
m OF (to) 
T eff dm 

where we introduced F = F — T e ffS c with 

F (to) = - lim — log / DyjZ m [$] . 
g^O to J 



It is now possible to integrate out the auxiliary variable 
ip, yielding an m-times replicated theory of the original 
variables <p with infinitesimal inter-replica coupling: 




similar to a random field model with infinitesimal ran- 
domness g. The major difference here is that in systems 
with a tendency towards self-generated glassiness, the ini- 
tial infinitesimal randomness g will self consistently be re- 
placed by an effective, interaction induced, self-generated 
randomness. Specifically this will be the off diagonal el- 
ement of the self energy in replica space. 

At this point it is useful to discuss similarities and 
differences of the present approach if compared to the 
conventional replica approach of systems with quenched 
disorder. First, on a technical level, the replica index has 
to be analytically continued to to = T/T e g < 1 and not 



to zero. This reflects the fact that slow metastable con- 
figurations do not equilibrate at the actual temperature, 
T, but at the effective temperature T eff . In other words, 
the system has an essentially equal probability to evolve 
into states which are spread over a spectrum with width 
T eff , even as T — > 0. If one applies the present approach 
to a system with explicit quenched disorder, where one 
can apply the conventional replica theory, and assumes 
replica symmetry, it turns out that to corresponds to the 
break point of a solution with one step replica symme- 
try breaking of the conventional replica approach and 
both techniques give identical results. Alternatively one 
might also consider the model, EqQJ in the limit of an in- 
finitesimal random field with width g and break point of 
a one step replica symmetry breaking to solution. This 
implies that the present approach captures the essence 
of glassy behavior in systems with two very distinct typ- 
ical time scales. This will become particularly clear if 
we compare our results with the one obtained within the 
solution of the dynamical Schwinger-Keldysh theory in 
the appendix. A more physical relationship between the 
two replica approaches can be obtained by realizing that 
the typical free energy of a frozen state, F, can also be 
written in the usual replica language via 

F = lim - ( - 

two n V ^ J 

where the average is performed with respect to the dis- 
tribution function p^ oc Z™ . The distribution function 
of the self generated randomness is non-Gaussian. For 
example if to — > 1, e~ 5 ^ v l can also be interpreted as the 
generating functional of the distribution p^w^ It is be- 
cause this distribution is characterized by "colored noise" 
that we find a self generated glassy state of the kind 
discussed here. Finally, if one replaces ([V 2 + Qq] <p) 
in EqQ] by the usual (V<^) 2 -term (which is the go — > 
limit after appropriate rescaling of <p, c, u. etc.) there 
is no glass asj-»0, making evident that self-generated 
glassiness is ultimately caused by the uniform frustration 
of the finite qo problem where modulated configurations, 
ip (x) oc cos (qo • x), with |qo| = qo and arbitrary direc- 
tion have low energy. 

IV. DYNAMICAL MEAN FIELD THEORY 

Due to the mean field character of the theory, it is 
appropriate to proceed by using the ideas of the DMFT 
for equilibrium many body systemsii^i^ and assum- 
ing that the self energy of our replicated field theory is 
momentum independent. Physically, ignoring the mo- 
mentum dependence of the self energy might be justified 
by the fact that a glass transition usually occurs in a 
situation of intermediate correlations, i.e. when the cor- 
relation length of the liquid state is slightly larger but 
comparable to the typical microscopic length scales in 
the Hamiltonian. 1 ^ The free energy of the system can 
be expressed in terms of the Matsubara Green's function 
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Figure 1: DMFT Phase diagram for the model Eq. 1 , within 
the one-loop self consistent screening approximation. The 
solid line represents Ta whereas the dashed line represents 
Tk. The inset shows the low temperature region of the phase 
diagram where Ta and Tk have merged and the effective 
temperature at the transition is larger than T (»ti(Ta) < 1). 
The line is simply a guide to the eye whereas the dots are 
the actual numerical result. The numerical simulations in 
the paper where carried out using the parameters ro = —6, 
go = 0.3 and u = 2ir 2 \r \. The stars refer to the (T, c) points 
of which the results in Fig. 4 are shown. 



Gff (ui n ) = ((p^(uj n )(p b _ cl (-0On)) and the corresponding 
local (momentum independent) self energy T, a/3 (ui n ) as 

F (to) = tr (EG) - tr log G + $ [G] , 

where the self energy is given by E = ^fg^ , and $ [G] is 
diagrammatically well defined for a given system^. 

The difference from the usual equilibrium DMFT 
approac h 17 i 18 i 19 is the occurrence of off diagonal elements 
in replica space, allowing us to map the system onto a 
local problem with the same interaction but dynamical 
"Weiss" field matrix. This brings us to an effective zero 
dimensional theory similar to the mode-coupling theory 
of classical glasses^ We furthermore make the ansataSS 

GfW=5 q k)^ + J q ^, (4) 

with static off diagonal elements. A similar ansatz for the 
self energy leads to the following two Dyson equations for 
the diagonal and off diagonal propagators: 



Gq (W n ) = [£ 0: q - S <? 1 > 

•Fq = ^Gq (0) [Sq 1 (0) - TO^/T] 



(5) 



Glassiness is associated with finite values of the Edwards- 
Anderson order parameter, J-q, whereas for J-q — we 
recover the traditional theory of quantum liquids. The 
structure of the Dyson equation already gives us crucial 
informations on the nature of the glass transition. From 



Eq0 follows immediately that, contrary to the classical 
case where glassiness can occur with to = l 1 ^' 14 , in the 
quantum limit (T — > 0) the only way Tq can be non-zero 
is to have to — > such that m/T = 1/T e g ^ 0. This 
imposes a constraint on the replica symmetry breaking 
structure in the quantum limit of the replica approach 
developed in Refl2fi Also, it is clear that defines a 
new length scale of the problem that is associated with 
the glass transition^ In the classical glass transition the 
relation m = 1 is satisfied at Ta- Then the two Dyson 
equations can be decoupled into an equilibrium, diago- 
nal (in replica space) part, and a non-equilibrium, off- 
diagonal part. This allows us to interpret Eg as related 
with the equilibrium correlation length and to associate 
Ejr with the Lindemann length, associated with the typ- 
ical length scale of wandering of defects of the equilib- 
rium structure, see Refi^ for details. However, it follows 
from the structure of the self-energies that, in the quan- 
tum limit, where m < 1 at Ta, the two Dyson equations 
cannot be decoupled anymore. In this case, the two self- 
energies will combined define a correlation length and a 
Lindemann length which are not independent, but rather 
closely intertwined. 

We solved the impurity problem within a self consistent 
large- TV approach, i.e. we generalize the scalar field <p to 
an iV-component vector and consider the limit of large-iV 
including first l/7V-corrections. This approach was used 
earlier to investigate self generated glassiness in the clas- 
sical limi t 1 ^ 14 . In this limit it is also possible to solve the 
DMFT-problem exactly2£, demonstrating that glassiness 
found in the approximate large-iV limit is very similar 
to the exact, finite N theory and thereby supporting the 
applicability of the large TV-expansion. 

Within the self consistent large- N approach, the diag- 
onal and off diagonal element of the self energy are given 
as 



Eg (oJ n ) 



DfQ («„)) 
+T^D g (w n + U) m ) Q (u m ) 

7) I 

2 



N 



(6) 



with Hartree contribution 



yH 



-uT ^2 G ( w « 



uJ- 



as well as 

Dg (ui r , 



Ug {u n ) + 2jFg {uo n ) 



D 



and 



T 



bubble 



-l^D 2 g (0) 

1 + jJ^Dg (0) /T eff 



(7) 



diagram 

Q 



i T T, m G{uj 

turn averaged propagators. 



J-0jGq(uJ n ) and T = j 



Ug (u n ) 

Here Q (w„ 



d 3 q 



rJ ^sJ 7 q are the momen- 
The DMFT is usually 



formulated on a lattice and it is possible to chose the 
same dimension (e.g. inverse energy) for the momentum 
dependent and momentum averaged propagator, by 
assuming the lattice spacing equal unity. In a continuum 
theory the role of the lattice spacing is played by the 
inverse upper cut off of the momentum integration A, 
which enters our theory in Eqs0 through the constant 
7 = A~ 3 . In the case of the Hamiltonian, Eq^ all 
momentum integrals are convergent as A — > oo and 
the scale which replaces the cut off is qo , leading to 
7 = qQ 3 . We solved this set of self consistent equations 
numerically. Before we present the results we must 
discuss the stability of the ansatz, EqQ] 

The local stability of the replica symmetric ansatz, 
EqQ] is determined by the lowest eigenvalue of the Hes- 
sian matrix SG auf G -,i, ■ Proceeding along the lines of 
Ref 36 and diagonalization over the replica indices leads 
to the following matrix in momentum space 

Mq.q, =<5(q-q')<? q 2 (0)+C, (8) 

where 

8 2 $ 6 2 $ 5 2 $ 

SGal 3 SGaf 3 ~ SG a ^SG aS + SG^SCr 5 ' ^ 

with distinct a, f3, S and S. Diagrammatically, C is a sum 
of diagrams with four external legs with at least two dis- 
tinct replica indices. Thus, as g — > the constant C 
vanishes if Tq — > 0, an observation which will be relevant 
in our discussion of the nature of the zero temperature 
glass transition as a function of c. We find for the lowest 
eigenvalue, A, of Mq,q>: 

r- 1 + *)-'■ m 

If A > the mean-field solution is stable, whereas it is 
marginal for A = 0. The fairly simple result, EqEl for 
the stability of the replica structure is a consequence of 
the momentum independence of the self energy within 
DMFT which guarantees that $ [G] only depends on the 
momentum averaged propagators. Otherwise C would 
depend on momentum, making the analysis of the eigen- 
values of Mq,q> in Eq|H|much more complicated. Thus, 
the investigation of our problem within dynamical mean 
field theory is not only convenient to obtain solutions for 
Gq (u> n ) and J-q, it is also crucial to make progress in the 
analysis of the stability of this solution. We expect that 
it will be impossible to find a stable replica symmetric 
solution, EqQ] once one goes beyond DMFT. 

In real physical systems the slow degrees of freedom 
relax on a finite time scale r Q and the effective tempera- 
ture T e ff = T/m depends not only on the external param- 
eters like T and pressure, but also on the cooling rate, 
or equivalently on the time, t w , elapsed after quenching 
and thus on the events which, for a finite range system, 
can take the system to different states of the spectrum 
characterized by T e g. Accounting for these effects goes 
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Figure 2: Low temperature behavior of the Edwards- 
Anderson order parameter as a function of the intensity of 
quantum fluctuations. The temperatures are in the region 
displayed in the inset of Fig. Q Inset: dependence of the 
order parameter at the transition temperature T\. 

beyond a mean field treatment. However, for t w <C r a 
mean field theory should apply. In fact, since within the 
mean field approach the time r Q is infinite, this is in- 
trinsically the regime we are constrained to within the 
mean field treatment. Then, the most important con- 
figurations of the order parameter are those which allow 
the system to explore the maximum number of ergodic 
regions. The best way to achieve this, without being un- 
stable, is through interconnecting saddle points^ This 
leads to the marginal stability condition, A = 0, which 
we use to determine the effective temperature T e ° ff , which 
corresponds to the effective temperature right after a fast 
quench into the glassy state, i.e, T e ° ff = T e g(t w <§; r a ): 

c{n*Y l = j^ii^x*)- (id 

Within the large- N approximation used to determine the 
propagator we can also evaluate the constant C of EqEJ 
leading to 

C{T° s ) = ^[2Dl(0)Ilr-Vr} (12) 

which is an implicit equation for T e ° ff . Note that, as we 
argued before, C vanishes for T — » 0. Moreover, since the 
integral on the right hand side of EqEJis bounded from 
above, C and therefore Tq must vanish discontinuously at 
the transition even for T = 0. As explained above, from 
the structure of the Dyson equation in replica space it 
immediately follows that for c larger than some value ca, 
T e ° ff /T must also jump discontinuously from T® s /T > 1 
in the glass state to T e ° ff /T = 1 in the quantum liquid 
state. 

Together with EqEJwe have a closed set of equations 
which describe the quantum glass within mean field the- 
ory. We have solved this set of equations for the model, 



Eq^ and indeed found that there is a glassy state be- 
low a critical value for c. Most interestingly, for increas- 
ing c (i.e. increasing quantum fluctuations) the rapidly 
quenched quantum glass, which at some point becomes 
unstable, undergoes a discontinuous reorganization of the 
density matrix upon entering the quantum liquid state. 
It is crucial to solve this many-body problem within some 
conserving approximation, i.e., based upon a given func- 
tional $(G), which must be simultaneously used to de- 
termine S and C. The set of equations was solved numeri- 
cally. The Matsubara frequency convolutions were calcu- 
lated on the imaginary time axis via Fast Fourier Trans- 
form algorithm using 2 15 frequencies. This accuracy is 
needed mostly to be able to find solutions of EnslTTland 

[H. 

The transition line between the liquid and glassy states 
in the (c, T) space is presented on Fig^ Note that the 
low temperature (quantum regime) behavior is qualita- 
tively different from the classical stripe glass. In the 
quantum limit, Ta and Ik merge and the effective tem- 
perature at the transition T e ° ff (Ta) is always larger than 
Ta, i.e., m < 1. Due to the reentrant character of the 
transition the quantum glass can also be reached by heat- 
ing up the system. By generalizing Brazovskii's theory 
of the fluctuation induced first order transition to the 
quantum case, we found a similar reentrance behavior for 
this equilibrium transition, suggesting that this peculiar 
shape of the phase border is determined by the increas- 
ing relevance of fluctuations with wave vector q = [c?o | 
as one crosses over from a quantum to a classical regime 
(see corresponding remarks made in sec. II). Another way 
the quantum glass can be reached is of course via a "c 
quench". Note also that the reentrant behavior we find 
within our approach happens at the point where numer- 
ically, at the same time, Ta — Tk vanishes (to within 
numerical precision), to (TO starts falling with a larger 
derivative (see the inset of FigEJl, and J- (T^) reaches 
a minimum (see inset of Fig. [2j before plateauing as 
T — > 0. This suggests that the reentrance behavior and 
the change of character of the transition are closely re- 
lated. In Fig. 3 we show the dependence of the replica 
symmetry index to = T/T e ° ff as function of c for different 
temperatures. 

In the classical limit, c — > 0, T<± changes discontinu- 
ously at Ti, whereas T e ff changes continuously (m = 1), 
i.e. the relevant metastable states - within which a clas- 
sical glassy system gets trapped into - connect gradually 
to the relevant states which contribute to the liquid state 
partition function. A similar behavior in the quantum 
limit would imply that T e g and therefore JF q and C van- 
ish continuously at the quantum glass transition. How- 
ever, as discussed above, this is not possible if one uses 
EqEJto determine the effective temperature. Thus, T e ff 
changes discontinuously and one might expect a nucle- 
ation of liquid droplets within the unstable glass state to 
be important excitations which cause a quantum-melting 
of the glass. Even though one can formally introduce, 
along the lines of Ref ill a latent heat SQ = T e ffS c at this 
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Figure 3: Dependence of the effective temperature as function 
of c for different temperatures. Inset: dependence of T® s /T\ 
on Ta in the low temperature regime displayed on the inset 
of Fig. [I] 

transition, we do not know of a scenario which, within 
mean field theory, allows this energy to be realized within 
the laboratory. The appearance of a first order like tran- 
sition as one enters the quantum regime, however without 
reentrance behavior, was first pointed out in Reflifll in a 
related case of spin glasses with quenched disorder. 

Finally, we make contact between our theory and the 
Schwinger-Keldysh approach used in Ref l2flL which gives 
a set of coupled equations for the symmetrized correla- 
tion function C x , x / (t, t w ) = \{[<p x (t w + r) , ip x > (t w )] + ) 
and the retarded response function G x x , (r, t w ) = 

-iQ(r) ([ip x (i w +r) (*w)]_), where [, ] ± = AB ± 
BA . In the classical limit C x . x ' (r, i w ) = 
(Vx (t w + t) ip x > (i w ))as usual. If t w , is comparable to r 
the dynamics is complex and depends on the nature of 
the initial state (aging regime). On the other hand, for 
t w large compared to r, C (r, i w ) and G r (r, i w ) are ex- 
pected to dependent only on r (stationary regime). Cor- 
respondingly, one can decompose the correlation function 
into aging and stationary contributions 

C(r,t w ) =Cag(t,M + Cst(t) , (13) 

and similarly for G r (r, f w ). Cugliandolo and Kurchan^i 
showed within the mean field theory of classical spin 
glasses that one cannot decouple the stationary dynamics 
from the aging regime. Instead, the system establishes a 
"weak long term memory" and one has to solve for the en- 
tire time dependence. An elegant way to encode the aging 
dynamics is a generalized fluctuation dissipation theorem 
(FDT) G AG (r, t w ) = -T-^Cag (r, t w ), with effective 
T e ff. This approach was generalized to the quantum case 
in RefE3 

As shown in the appendix, we find a complete equiv- 
alence of the Schwinger-Keldysh theory of Refl2fil. if 
applied to the quantized Brazovskii model, with our 
replica approach if we identify (after analytical con- 
tinuation to real time): Q (t) — Gg T (£)and T = 
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lim tw _ 00 lim^oo Cag (t w ,t). T e ° ff , which follows from 
marginality, EqEl is identical to the one of the gen- 
eralized FDT, when we assume t w <^r a £& 



V. ASPECTS BEYOND THE STRICT MEAN 
FIELD THEORY 

In this section we discuss several physically significant 
conclusions one can draw from our theory which go be- 
yond the strict mean field limit. In particular we will 
discuss several aspects related to dynamical heterogene- 
ity in glasses. 

Going beyond mean field theory, for t w 3> r Q , the 
marginal stability cannot be sustained, since the system 
can save free-energy via droplet formation^ which drives 
the system towards equilibrium. As discussed in RefEfj 
the free-energy gain is due to a gain in configurational en- 
tropy which is inhibited within mean field theory due to 
infinitely large barriers, but possible in finite subsystems 
where transitions between distinct metastable states be- 
comes allowed. Thus, the glass might be considered as 
consisting of a mosaic pattern built of distinct mean field 
metastable states. The size of the various droplets form- 
ing the mosaic is determined by a balance of the entropic 
driving force (proportional to the volume of the drop) 
and the surface tension^. 

One way to account phenomenologically for such be- 
havior within our theory is to assume that T e g becomes a 
time dependent quantity^ and the exploration of phase 
space allows the system to "cool down" its frozen degrees 
of freedom by realizing configurational entropy. One 
would naturally expect T e g(t w ) to decrease towards a 
value T°g until either limbec S C {T^)/N = or T£§ = 
T. In the former case there is no extensive entropic driv- 
ing force anymore which favors the exploration of phase 
space, whereas in the latter case the system has reached 
equilibrium but with, in general, finite remaining configu- 
rational entropy. These two regimes are separated by the 
Kauzmann temperature where limjv^oo S C (T^)/N = 
at = T simultaneously^. We cannot, of course, calcu- 
late the explicit time dependence of T e g(t w ) here. How- 
ever, we can parametrically study S c (T e g) versus T e g, i.e. 
keep the replica value m an open parameter of the the- 
ory and analyze whether the trends for the variation of 
S c (T e g) at fixed temperature T are sensible. Since we de- 
termined m previously by the marginality condition, an 
effective temperature T e g < T e ° ff implies that the replicon 
eigenvalue is different from zero. We found numerically 
that A < 0, reflecting an instability of our solution likely 
related to some kind of dynamical heterogeneity. We ar- 
gue that this heterogeneity is different for temperatures 
close to Ta and Tk- 

We first discuss the behavior below but close to Ta in 
the classical regime where marginality gives a continuous 
change of m, i.e. (Ta) = T. In this regime the sys- 
tem is close to equilibrium and it should be sensible to 
consider small fluctuations around the mean field solu- 




1 1.1 1.2 1.3 

T e f / T 

Figure 4: Parametric plot of the relaxation of T e s and S c 
for varying t w after quenches to the 3 points indicated by 
stars in the phase diagram Fig[2 In addition S c is shown for 
one point right at the quantum glass transition with Ta = 
0.01 > T. The direction of evolution of the waiting time t w is 
indicated by the arrows. The system can relax from a state 
with marginal stability with T e ° ff to 3 different final states as 
explained in the text. Inset: configurational entropy at the 
transition point as function of Ta. 

tion. Such small fluctuations are then dominated by the 
eigenvectors of the replicon problem which correspond 
to the lowest eigenvalue. From EqOD we can easily de- 
termine the momentum dependence of this eigenvector, 
which can be interpreted as the fluctuating mean square 
of the long time correlation function J-q away from its 
mean field value. It is given by 

* q = *o-^- (14) 

where, Wo is a normalization constant. At marginality, 
these modes in correlation function space are massless 
and thus easy to be excited. The typical length scale 
of these correlation function fluctuations are determined 
by Gq=q , i.e. are confined to a length scale determined 

by Tig 1 / 2 . Since this length is not the actual correla- 
tion length, but rather determined by the shorter Linde- 
mann length discussed on Refill the wandering of de- 
fects seem to set the scale on which the dynamics close to 
Ta evolves. Thus correlations over the Lindemann length 
are fluctuating in space. Even though the fluctuating 
object is characterized by a rather short scale, its fluc- 
tuations are, of course, correlated over larger distances, 
as characterized by the nonlinear susceptibility ^ a/3 ^ s = 
((p a (x) tfif3 (x) ifiy (x) ips (x))- x a ^ 1& is the inverse of the 
Hessian and thus diverges if A — > 0. Thus we conclude 
that close to Ta this heterogeneity is driven by the cor- 
relation function fluctuations of shape Wq of Ea ll4l We 
believe that this "Goldstone" - type heterogeneity is very 
similar in character to the recent interesting approach to 
heterogeneity resulting from the assumption of a local 
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time reparametrization invariance 3 ^. This is supported 
by the close relation between the reparametrization in- 
variance and marginality as shown in Reflifl In addition, 
the approach of Reflsfll considers fluctuations relative to 
the mean field solution with stiffness proportional to the 
(mean field) Edwards- Anderson parameter. Thus, small 
fluctuations caused by the marginality of the mean field 
solution are considered, similar to the ones given in Eall4l 
Further away from Ta such a linearized theory is likely 
to break down because additional eigenvectors, not re- 
lated to the marginal eigenvalue, become relevant and 
non-Gaussian fluctuations come into play. This should 
always be the case where m is not close to unity, i.e. 
in the classical regime for temperatures below Ta and 
everywhere in the quantum regime. One might expect 
droplet-physics to become important ther*2IL. 

In both cases it is useful to analyze the evolution of 
the spectrum of states of this formally unstable theory, 
particularly if the replica structure of the theory remains 
unchanged. As shown in Fig.l, we found numerically for 
the model, EqQJ that for t w » r Q , three different possible 
final situations, result, depending on the relation between 
the bath temperature T and the Kauzmann temperature 
Tk'- (1) Ta > T > Tk, T e ff will relax until it reaches 
T, but an excess entropy will remain. (2) T = Tk, T e ff 
will relax until it reaches T, in a state with zero con- 
figurational entropy. (3) T < Tk, T e g will relax until 
all the excess entropy vanishes, but the system remains 
in a non-equilibrium state with > T. In this case 
there are still many (even though less than exponentially 
many) states, distributed according to a Boltzmann func- 
tion with an effective temperature, T°jf. Their energies, 
fy, differ by non-extensive amounts in a range of order 
T°g < T e ° ff . At a critical value of the quantum param- 
eter slightly below ca we find that Ta and Tk merge 
and the nature of the glass transition changes. The sys- 
tem is either in a quantum fluid or in a non-equilibrium 
frozen state and the identification of the glass transition 
using equilibrium techniques alone becomes impossible. 
Note however that in the quantum glass regime, even at 
temperatures T arbitrarily close to Ta, for t w 3> t q we 
always obtain Um^ w S C {T™)/N = for T™ > T (see 
dashed curve on Fig. 0}. 



VI. CONCLUSION 

In summary, we have presented an approach to self 
generated quantum glasses which enables the counting 
of competing ground state energies (quasi-classically of 
long lived metastable states) in interacting quantum sys- 
tems. Technically very similar in form to the traditional 
quantum many body theory of equilibrium systems, it 
allows one to investigate whether a given system exhibits 
self generated glassiness as a consequence of the frustrat- 
ing interactions. Slow degrees of freedom are assumed 
to behave classically and are shown to equilibrate to an 
effective temperature which is nonzero even as T — > 



and which characterizes the width and rigidity of the en- 
ergy landscape of the competing states of the system. 
Applied to the specific model, Eq^ we do find a glass 
below a critical value for the quantum fluctuations. Using 
a marginality criterion to determine T e g, we can gener- 
ally show that quantum glass transitions are bound to 
be discontinuous transitions from pure to mixed quan- 
tum states. This leads to the interesting question of how 
quantum melting of the non-equilibrium quenched states 
can occur via nucleation of the corresponding quantum 
liquid state. Even going beyond mean field theory by as- 
suming a time dependent effective temperature, we find 
that T e g saturates (at least for extremely long times) at 
a value T°g > T. Finally we made connection to the 
dynamical approach for non-equilibrium quantum many 
body systems of R,ef l2flL which shows that our theory 
properly takes into account the effects of aging and long 
term memory. We believe that the comparable simplic- 
ity of our approach allows it to apply our technique to a 
wide range of interesting problems in strongly interacting 
quantum systems. 
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Appendix A: SCHWINGER-KELDYSH 
FORMALISM 

In this section we apply the Schwinger and Keldysh 
close-time path Green function formalism applied to 
quantum glasses without quenched disorder. We follow 
closely Cugliandolo and Lozano2^, who applied the tech- 
nique to spin glasses with quenched disorder. In this for- 
malism response and correlation functions are treated as 
independent objects, coupled through a set of equations 
called Schwinger-Dyson equations. A perturbation the- 
ory scheme can be set up by considering the generalized 
matrix Green function 

g ,(tt')-( G x,x<M'n (A1) 

and self-energy 

g , (tt >)-( ^ X ,( M ')\ 
x - x 1 ' j ~ \ ^(t,t') E x , x ,(i,f') ) ■ 

where 

Cx,x'(i,t') = ^([<Px(tw+T),(£> x '(£w)] + ) 

<")(*,*') - T i0(±(i-*'))([M*),¥V (*')]> 
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In analogy to conventional perturbation theory, the 
components of G and E obey the Schwinger-Dyson equa- 
tions 



when accumulated over long times, gives a finite contri- 
bution to the dynamics. Thus, following along the lines 
of Ref S3, we get for the for the aging regime 



G q = G' 0q + G5 q ® E q ® G q (A2) 
Gq = Gq ® [G£ q 1 ® G 0q ® G° q 1 + E q ] ® G° .(A3) 

We use ® to distinguish the matrix product over 
time (the time convolution) where A <E> B (t, t') 
j dsA(t, s)B (s,t') from the scalar (element-wise) prod- 
uct, where AB(t,t') =A(t, t')B(t, t'). 

The first Schwinger-Dyson equation IA2I is similar to 
conventional perturbation theory. However it is coupled 
to the second equation of C {ti^t^j which admits non- 
trivial solutions. For simplicity, we choose an initial con- 
dition such that Go = and ImGg = 0. This leads to 
thermalization of G r (fulfillment of the FDT) in the ab- 
sence of non-linearities, so it is an appropriate condition. 
Let us call £2 = t w and t\ = t w + r, i.e., the correla- 
tions between instants separated by r are measured after 
the waiting time t w has elapsed. The glassy dynamics ap- 
pears in a regime of t w — > 00 and r — ► 00 . To proceed we 
assume that in this limit all correlation functions can be 
decomposed into a slow, non-time translation invariant, 
aging part (AG) and a fast, time translation invariant, 
stationary part (ST) like in Eall3l 

In this approach the stationary term decays on a char- 
acteristic time scale Tp and it represents the correlations 
between degrees of freedom which are in equilibrium with 
the thermal reservoir at temperature T, i.e, retarded and 
Keldysh correlation functions are related by the FDT. 
The aging part on the other hand depends weakly on t w 
(aging phenomena) and varies slowly on r in a character- 
istic large time scale r Q 3> 173 , which allows us to neglect 
its t derivatives. Moreover, as in reference HE we en- 
force a relationship between G r AG and Gag by defining 
an effective temperature T e ff at which the long time cor- 
relations thermalize through a generalized FDT relation 



Gag q ( T ; *™ ) — ^ ag ( T ) t w ) <?st q (0) 



(A6) 



Analogously, the second Schwinger-Dyson equation 
l!A3j) can be solved for the same two regimes, yielding 



GfTqM 



and 



Gag q (t, t v 



^STM^STqM" 



Sag (t, tw) Ga,Tt 



^STq (°) lfe S AG ( T ' *«») 



(A7) 



(A8) 



Note that if we identify the real time correlation functions 
with the inter-replica correlation functions T and Q as 
follows 



(A9) 

u + i5) (A10) 



lim lim GagCM™) = T 

T — >00 t w — >00 

0ST M = Q 



we get that the Schwinger-Dyson equations IA5I and 
IA8I are exactly equivalent to the replica equations in 
provided that the self energies in the two schemes 
are also equivalent, i.e., Sag = — and Eg T (u>) = 
Hg (iu> n —>■ lu + iS). To prove this last requirement, we 
apply the same one-loop self-consistent screening pertur- 
bative scheme for the real time DMFT self energies E, 
This gives 



E = 





D r + D> 



C®D r + D®G 
where dressed interactions are given by 



C ® D a + D <g> G a 
D a ®G a + D r ®G r + D®C 



D r = u + vH r ®D r 
D = uIT ® D + uli c 



D a , 



(All) 
(A12) 



Gag (r,t) 



d 



T eS (t) dt 



Gag (r,t) 



(A4) 



For t <Si r a we can neglect the time dependence of T e g, 
leading to a constant T e ff (t) — > T e ff. 

Using the above assumptions, the first Schwinger- 
Dyson equation ilA2j) in the stationary regime can be 
solved by a Fourier transform, which gives 



^STq M ~ r r-l 7 a 1 



(A5) 



Here we will explore solutions where G AG (r, t w ) is small 
but finite. This corresponds to the weak long term sce- 
nario where, even though lim T „ +00 lim^oo G AG (r, t) — > 
0, the integral lim^oo J dt'G r AG (r,f) is still finite. In 
other words, we consider that the system keeps a vanish- 
ingly small memory of what happened in the past which, 



and the polarization bubbles are given by the element- 
wise products (no time or space convolutions) 



IT 



7 GG r 



n = 7 



(G a ) + (G r y + G 2 



(A13) 
(A14) 



The stationary and aging contributions to the self ener- 
gies and polarizations can be calculated analogously by 
using the definitions l|A13Jl and l|A14l) and separating the 
stationary from the aging contribution according to their 
asymptotic time behavior. This gives 



yr 
^ST 



Sag — 



Gag ® £>st + -Dag € 

+ G ST ® £>ST + D ST 

Dag ® Gag 



Gr 
ST 

$ Gg T 



and 



n 



ST 



7 (GAcGg T + GsTGg T ) 



(A15) 



(A16) 



(A17) 
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and 



IIag = 7Cag 



(A18) 



It is easy to verify accordingly that n ST and EIg T are re- 
lated by FDT at the temperature T and IIag and II AG 
are related by the classical FDT at a temperature T e ff • 
The last term in En I A 151 is exactly the Matsubara con- 
volution on 03 . The remaining terms are calculated as 
follows: 



v r ST H 



+ n ST (w) 



D 



AG 



-n AG v r ST (lo = o) 
v r ST ( w = o)" 1 + ^n AG 



(A19) 



(A20) 



Using the relations IA9I and IA10I and performing 
the usual Matsubara sums we get,HV = IIag and 
Ilg (iui n — ► u) + iO) = IIg T (u>) which yields Dag = Dyr 
and-Dg {iui n — > cj + iO) = Pg T (w). Therefore, equa- 
tions ijAlflj) . (|A2np> . HAX.qfl . and <|A16j) together prove 
that Eag = Ejf and E ST (ui) = Eg (iw n — > u + iS)and 



that there is a complete connection between replica 
and Schwinger-Keldysh formalisms. There is still one 
more independent equation, namely, Eq. lAfU in the 
Schwinger-Keldysh formalism which bares no analog on 
the Dyson equations for the inter-replica correlation func- 
tions. However, if we integrate over q both sides of lAfil 



substitute the definition E 



AG 



C AG ®D r AG +D AG ®G 



into El with D r AG = II AG X>| T (0)we obtain 



AG 



[22>§ T (o)n£ G -D£ G ] 



9 C7 S t'(0) = 1 (A21) 



(27T) 



which, is exactly the marginality condition of the replica 
approach. This proves our previous statement that the 
saddle point condition gives the dynamical behavior at 
the time scales t w — > oo, r — > oo with t w <C r Q . In 
this limit the effective temperature is not yet in equi- 
librium. In the opposite limit one cannot disregard the 
time dependence of T e g- (i)and it is thus not possible to 
write a closed form like equation l!A21j) . Nevertheless, 
the equilibrium approach based on the replica trick, to- 
gether with the droplet relaxation picture enables us to 
access the behavior of T e ff (^w) m the time limit t w 3> t q . 
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